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Abstract 

A model for the 7p — ► Tr + ir~p reaction developed earlier is extended to 
account for all isospin channels. The model includes N, A(1232), iV*(1440) 
and iV*(1520) as intermediate baryonic states and the p-meson as an inter- 
mediate 2ir resonance. Although many terms contribute to the cross section, 
some channels exhibit particular sensitivity to certain mechanisms of reso- 
nance excitation or decay and the reactions provide novel information on such 
mechanisms. In particular the jN — > N* (1520) — > An process affects all the 
channels and is a key ingredient in the interpretation of the data. Comparison 
is made with all available data and the agreement is good in some channels. 
The remaining discrepancies in some other channels are discussed. 
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1 Introduction. 



There are three possible double pion photoproduction reactions on the proton, and 
three on the neutron: 

(a) jp — ► 7T + 7T~p (d) 7n — > 7T + 7T~n 

(b) jp — > 7T + 7T ?T, (e) '-fU —>■ 7l~7T°p (1) 

(c) 7J0 — > 7r°7r°p (/) 772 — > 7r 7r°n 

These reactions have been extensively studied experimentally in the past ([Jl], @, 
0. U> §, @> 0J). New improvements in experimental techniques and facilities have 
reopened the study of these reactions at Mainz, with two experiments on the proton 

SEP- 

From the theoretical point of view, only the reaction (a) has been studied with 
one early model JTC| ] which considers only 5 Feynman diagrams, and a more com- 
plete one [T^] which includes N, A(1232), iV*(1440) and iV*(1520) as inter- 
mediate baryonic states and the p-meson as intermediate 2n -resonance. This model 
reproduces fairly well the experimental cross sections below E 7 = 800 MeV, and 
the invariant mass distributions even at higher energies [|PT|]. 

On the other hand, the double pion photoproduction reactions on nucleons has 
been recently studied at threshold from the point of view of chiral perturbation 
theory [ |T3| , [RJ, [13], [R]]. In ref. [[RJ the authors show that due to finite chiral loops 



the cross section at threshold for final states with two neutral pions is considerably 
enhanced. 

Thus, with these reactions becoming a target of new experimental and theoreti- 
cal studies, plus interesting medium effects predicted for the (7, reaction in 
nuclei [|T7p (in analogy to those already found for the (71", 2n) reaction [|TSJ, [HJ), a 



thorough theoretical study of the 7iV — > -rniN reaction is timely and opportune. 

In this work, we have improved the model of reference [|TT]], and have extended 
it to the other isospin channels of Eq. (JT]). Moreover, although is not the purpose of 
this paper to perform a thorough analysis of the 7iV — > tttcN reaction at threshold, 
we have calculated the total cross section for the 7^ — > n + n~p and ^p — ► 7r°7r°p 
reactions at threshold energies, and compared our result with the result given in ref. 
| jT3| , [T3p based on Chiral Perturbation Theory. 

The model is based on the coupling of photons and pions to nucleons and reso- 
nances using effective Lagrangians and thus leading to a set of Feynman diagrams 
at the tree level. As well as in ref. [pj|, we do not implement unitarity in the final 
states, but we made an estimate of possible uncertainties owing to unitary correc- 
tions for the 7p — > n + ix~p channel. 



2 The model for the 7^V —> nnN reaction. 



2.1 the 7J9 — > 7r + 7r p channel. 



For the 7^ — > n + n p isospin channel we follow strictly the model of ref. [|TT|], 



with some slight modifications in the iV*(1440) and iV*(1520) couplings in order 
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to reproduce the new branching ratios of the last Review of Particle Properties [ ^0| ] 
(see Appendix A of the present paper for Lagrangians and coupling constants), and 
we have also improved the iV* (1440) propagator, modifying the energy dependence 
of the width to account explicitly for the iV*(1440) decay into one and two pions. 

We classify our diagrams in one point, two point and three point diagrams, ac- 
cording to the number of vertices in the baryonic lines (see Fig. 1). Our basic 
components are pions, nucleons and nucleonic resonances: For the baryonic com- 
ponents we consider N, A(1232, J n = 3/2+, I = 3/2), iV*(1440, J n = 1/2+ I = 
1/2) and iV*(1520, J n = 3/2", / = 1/2). The iV*(1440) played an important role 
in the nN — ► nnN reaction [ [18| , [19| ] in some isospin channels and for this rea- 
son, together with the fact that the iV* (1440) excited in the 7iV collision is on-shell 
for E 1 ~ 600 MeV , we pay a special attention to this resonance. The iV*(1520) 
has a particularly large coupling to the photons and proves to be an important in- 
gredient, mostly due to its interference with the dominant term of the process, the 
7iV — > Arc transition through the gauge Kroll-Ruderman term. Higher resonances 
have a weaker coupling to photons and do not interfere with the dominant term, 
hence their contribution is small, at least for photon energies below 800 MeV, 
Mainz energies, where our model is supposed to work. Because of the important 
coupling of the p-meson to the two pion system and the •yir system we have also 
considered terms involving the p-meson. These terms are only relevant at high 
energies but show up clearly in the two pion invariant mass distributions at these 
energies [jTT|]. 

With these considerations, the basic diagrams which we consider have the struc- 
ture shown in Fig. 1. In diagram (a) and (b) the one point NNirn coupling stands 
for the s-wave nN interaction. We consider there only the isoscalar part of the 
amplitude. The isovector part is mediated by p exchange [^TJ] and hence it is ex- 
plicitly taken into account in diagrams (f) and (h). Diagram (c) contains the gauge 
term NNtt'j or Kroll-Ruderman term. We use a pseudovector coupling for the 
NNn vertex and this allows us to consider exclusively positive energy intermediate 
states in the hadronic propagators []22|]. In the two point and three point diagrams 
we include nucleon and the resonances as intermediate states. However, while all 
possible diagrams with N and A(1232) intermediate states are considered, we omit 
some with iV*(1440) intermediate states which are very small. For the iV*(1520) 
intermediate states we keep only the term which interferes with the dominant term 
of the amplitude (A(1232) Kroll-Ruderman term). In addition, all different time 
orderings of the diagrams are considered. 

The diagram (g) involves a gauge term pmcy coming from minimal coupling in 
the pun vertex which contains a derivative coupling. The diagram (i) involves the 
7p7r vertex which appears in the p — > 77? decay. Finally diagram (j) contains the 
anomalous coupling 737r P3[|. 

The Feynman diagrams considered and detailed calculations can be found in 
refs. [§ §. 

In spite of the small contribution of the iV*(1440) terms in the reaction 7^ — > 
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n + ix~p (less than 2% [[TTj])[], in order to be consistent we have revised the value of 
the coupling constants involving the iV*(1440) so that one may fit the new branch- 
ing ratios of ref. [ p0| ] (see Appendix A of the present paper for coupling constants). 
This updating gives us very small differences compared to the values given in ref. 

& 

The main change appears in the iV* Ait coupling constant due to the fact that in 
the study of the decay of iV*(1440) into Avr of ref. [[IT]] we considered the A(1232) 
as a stable particle. Now we modify this in order to take into account the finite 
width of the A(1232). The fact that the A(1232) width is not small compared to 
the mass difference between the iV*(1440) and the A(1232) makes this correction 
advisable. 

Making use of the Lagrangian (A. 7), the decay width in the cm. system of the 
iV*(1440) into An, considering the A(1232) as a stable particle (zero width), is 
given by: 



N* 



9 an** [ * m A 1 4 



d p — , -p 5(m N * - E A (p) - Lo{p)) 



* (2n) 2 fi 2 J r E A (p)2uj(p)3 1 

(2) 

where [L, m A and m N * are the pion, A(1232) and iV*(1440) masses respectively; 



p is the momentum of the outgoing pion; E A (p) = \Jm A + p 2 and to{p) = 
\J [i 1 + p 2 are the A(1232) and pion energies. 

In order to account for the finite A(1232) width, we have to replace the delta 
function of energy conservation in Eq. (g) by 

1 j m : (3) 

7T m N * - E A (p) -w(p) + |r A (v/sZ) 



where T^y/sX) is the delta decay width, and ^/sX is the A invariant mass, y/sX 



Pa" 2 ~ Pa 2 (Pa = (p° a ,Pa) is the A four momentum). With this modification, 
and using for T A the expression (B.32) of ref. [|TTJ], we get g A N*w = 2.07 after 
fitting the iV*(1440) decay width into A7r to the average experimental value (87.5 

On the other hand, in the iV*(1440) propagator 
1 rriN* 



y/s^-m N * + %T N *(y/s^) E N *(p N *) 



(4) 

where y/s N * = \/p%* 2 — Pn* 2 is the iV*(1440) invariant mass (p N * = (p° N *,p N *) 
is the iV*(1440) four momentum), and IV {y/s^) is the iV*(1440) decay width. 

We have included in the iV*(1440) decay width, both the decay into Nir and 
Ntttt, thus 



N" 



W $N*) — ^N*^Nn (V $N*) + T N *^ N7T7r (y/SN*) (5) 



1 However, for the case of photo-production of two neutral pions, the contribution of these terms 
is more relevant (between 10% and 60%, depending on the energy), thus it is important to treat them 
carefully. 
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where T N *^ N7T (^/s N *) is the A*(1440) decay width into Nn which is given by 

3 ( f\ 2 m 
2n \fi J m N * 

where / is the A*(1440)A7r coupling constant (see Appendix A), and \q \ is the 
pion momentum in the A* (1440) centre of mass system. We take a fraction of 
65% for the decay into Ayr [JZOfl. T N ^ Nnw (y/sJ^) is the A*(1440) decay width 
into Nitit, which we approximate by a constant, multiplied by three body phase 
space. This constant is fitted in order to reproduce the experimental iV*(1440) 
decay width into Nnn. Although part of the A* — > Nnn decay goes through the 
An channel (and this is explicitly considered in the Feynman diagrams) we follow 
the above prescription solely for the purpose of providing the energy dependence 
of the A* (1440) width. 

Apart from these small changes in the A* (1440) terms, the main improvement 
of the model of ref. [[IT]] is coming from the modification of the A*(1520)A7r 
coupling. In ref. [ |TT| ] we took the simplest Lagrangian compatible with the conser- 
vation laws: 

^att = ifNi*A^N>*<p X T x V A + h.c. (7) 

where ^ N >*, <P X and ^ A stand for the A*(1520), pion and A(1232) fields (A'* 
stands for the A* (1520) in the formulae from now on); T A is the transition isospin 
operator from 1/2 to 3/2 with the normalization in terms of a Clebsch-Gordan 
coefficient 

,3 w mt 1 . „ /l 3 

and /a^Att is the coupling constant, that we fixed from the data of A* (1520) decay 
into An of ref. [^4]. The sign however was chosen such as to have constructive 
interference below the A* (1520) pole. With the chosen sign the agreement with the 
data was relatively good, while with the opposite sign the discrepancies were about 
a factor of two and the qualitative features of the experiment were not reproduced 



-,m) = Cl-,l,-;m,u,M) (8) 



[11] 



Here, it is important to note that the A* (1520) can decay into An both in s- 
wave and d-wave [j2"0|]. However, the Lagrangian of Eq. (0) gives only the decay 
into s-wave, not d-waveQ. Then, we have modified the Lagrangian of Eq. ([7[) in 
order to get a term which gives a contribution to the <i-wave decay: 

Cn>*a* = Mn>* (fN»A* - 9 ~ E ^S\d l SjdA A T A ^ A + h.c. (9) 



2 In the previous version of the Review of Particle Properties [^4||, there were not branching ratios 
to s-wave and d-wave, and only the total decay into An was given. Thus, in ref. JIT] ] we assumed 
that all the decay was through the s-wave, and we neglected the d-wave. As the fraction of decay 
into the d-wave is even bigger than into the s-wave p0[|, corrections to the coupling of Eq. (f7|) must 
be now implemented. The influence of these changes in the results for the 7^ — * nnN reaction, 
although not large, certainly help to improve agreement with the data, as we shall see. 
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which gives us the vertex contribution: 

- i5H N » Av = - (f N ,* A „ + 9 Jt^& . $s ■ <f) T A (10) 

where S is the 1/2 to 3/2 spin transition operators with the same normalization as 
Eq. (|[) , n is the pion mass and q is thepion momentum in the iV*(1520) rest frame. 
In order to fit the coupling constants to the experimental iV*(1520) decay am- 



plitudes to An in s- and <i-wave PQ, [25Q, we make a partial wave expansion Q2q ] of 
the transition amplitude from a state of spin 3/2 and third component M, to a state 
of spin 3/2 and third component M': 

-i {^M'\8H N ^\\,M) = A S Y Q M - M \e,<j>) + 



A d C (2,~,~;M-M',M',M) Y 2 m - m '(9,0) (11) 



3 3 
2'2' 

where Y l m (9, <p) are the spherical harmonics, and A s and A d are the s- and d-wave 
partial amplitudes for the iV*(1520) decay into An. In order to compare with the 



experimental signs of the amplitudes p0| , |25| ] in Eq. < JllD we have followed the 
standard "baryon-first" phase convention []27|]. In this convention, spin couplings in 
all angular momentum Clebsch-Gordan coefficients are ordered such that the orbital 
angular momentum L comes before the intrinsic spin S. The angles in the spherical 
harmonics are those of the pion referred to the baryon. Then, making use of the 
coupling of Eq. ([T0|), the A s and A d amplitudes are given by: 

(12) 



A s = -a/47T [f N >*ATT + l9N'*Airf? 



A _ V^tt ~ q 2 
Ad — 3 9N'*Att ~ 

By using Eq. (01]) we obtain the iV*(1520) width into Att: 

IWa. = ^ — q (\A S \ 2 + \A d \ 2 ) (13) 

where m A and m N i* are the A and iV*(1520) masses respectively, and q is the mo- 
mentum of the outgoing pion in the N* (1520) rest frame. Now, we fit the iV* (1520) 



partial decay width in s- and <i-wave to the experiment [ pOQ (we take the average 
values for these partial decays, 8.5% for the s-wave, and 12% for the <i-wave). In 
addition, we impose the ratio A s /A d to be positive, as deduced from the experimen- 
tal analysis of the nN — > nnN reaction [|2"5|]. Thus, with this last restriction, we get 
two different solutions for the coupling constants which differ only in a global sign: 

(a) /jv* Att = 0.911 gN>*An = -0.552 

(b) f N '*An = -0.911 9N'*An = 0.552 

Now, the 7p — > n + n~p reaction allows us to distinguish between both solutions. 
In Fig. 2 we have plotted the total cross section for both solutions. As we can 
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see, only solution (a) fits the experiment, while the other one under-estimates the 
experimental cross section by a large amount. 

This strong dependence of the total cross section on the iV*(1520) coupling was 



already remarked in ref. [|TT|], and it is due to the important interference between 
the A Kroll-Ruderman term (diagram (z) of Fig. 3) and the •yp — > iV*(1520) — > 
An term (diagram (p) of Fig. 3) in s-wave [|TT|]. The amplitudes for these two 
diagrams have exactly the same structure, except by the iV*(1520) propagator and 
coupling constants (see Eq. (6) of ref. []1 where /n'^a-k would be substituted by 
JV'*Att + ^gwAn-q 2 / fJ> 2 )- Then, the iV*(1520) interference term is proportional to 
the real part of the iV*(1520) propagator jJlT]]. There is also a contribution to the 
total cross section coming from the imaginary part of the iV* (1520) propagator, and 
also from the Ad amplitude, which do not interfere. Thus, this interference is zero 
at ^Js = mjv*, where the real part of the iV*(1520) propagator vanishes. Then, 
what matters in the interference is the value of A s for different values of q than the 
one from the iV*(1520) on-shell. Then, the yp — > n + n~p reaction provides us with 
novel information about the q dependence of the A s amplitude, respect to the one 
obtained from the analysis of the nN — ► nnN reaction [[25|]. 

This q dependence in the s-wave amplitude is given by the chosen Lagrangian 
(on the other hand for the Ad amplitude there is no choice in the momentum de- 
pendence exhibited in Eq. (|12|)). It is, thus, possible to postulate other Lagrangians 
which would produce a different q dependence of the A s amplitude. For instance, 
in a preliminary work reported in [ |28"| , a different Lagrangian was investigated. 
There we proposed a similar Lagrangian to Eq. (Q), but with the 3/2 spin operators 
instead of the S l transition operators. With such a Lagrangian and an appropriate 
choice of parameters we were able to reproduce the experimental data, though the 
ratio A s /A d was negative. 

In order to investigate the most general q dependence of the amplitude, we re- 
place /jv'*A7r by 

jW* (l + e q2 ~ q y^ (15) 

where q is the momentum of the decay pion, and q 2 on _ she u is de momentum of the 
pion for an on-shell iV*(1520) decaying into An (\q m -sheii\ = 228 Me V). Then, we 
change e and compare the results to the data. We find that, to a good approximation, 
e = gives the best agreement with the data, hence supporting the Lagrangian of 
Eq.©. 

It is interesting to remark here that this q dependence of the amplitudes coin- 
cides exactly with the predictions of the non-relativistic constituent quark models, 
although the absolute values of the amplitudes are only qualitatively given by these 
models [p 



2.2 The others isospin channels. 

For the others isospin channels we use the same model as for the yp — > n + n~p case, 
changing the isospin factors and introducing some terms which are only relevant in 
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the case of neutral pions (although their influence in the cross sections is small). 
However, we have neglected some Feynman diagrams which were found to be very 



small in QTTQ. The Feynman diagrams that we take into account are shown in Fig. 
3, where all possible charge combinations allowed by charge conservation have to 
be implemented for each isospin channel. In the case of neutral pions and neutral 
deltas some of the diagrams in Fig. 3 are automatically zero. For instance, diagrams 
(c) and (d) of Fig. 3 with the photon attached to a neutral pion are zero. 

Nevertheless, some remarks are necessary for the isospin channels with neutral 
pions: When in the diagrams of Fig. 4 (a) and (b) we evaluate the contribution 
from the intermediate states of negative energy we get an extra contribution with 
the same structure as the Kroll-Ruderman term In the case of charged pions, 



we have neglected this contribution due to the fact that it is very small compared 
with the Kroll-Ruderman terms. However, for the neutral pions, as there is no Kroll- 
Ruderman term, this contribution has to be considered and, for a nucleon in the final 
state is given by 



/ 

- iT a = e — C(a) a ■ e (16) 

where 

C{a) = < (17) 
[ for neutron 

where (p°,p), {p'° ,p'), (q°, q) and (k°, k) are the four-momentum of the incoming 
nucleon, the outgoing nucleon, the pion and the photon respectively. 
We can write Eq. ( fRi| ) as: 

- iT = ~(1 + (N\r 3 \N))e f -Ca-e (18) 

where 

C=-q°( =r- + =r) (19) 

\p° + k° + E(p + k ) p' -k° + E(p'-k )J 

and r 3 is the third component of the isospin Pauli matrices: 

(pNp) = l ; (n\r 3 \n) = -l (20) 

The advantage of expression dl8[ ) is that it allows us to generalize it to the case 
where we have a iV* (1440) or a A(1232) in the final state by the following method: 
for the iV*(1440) we replace the coupling constant iViV7r (/) by the coupling con- 
stant iV*iV7r (/): 

- iT = i(l + (N\t 3 \N)) e f -Ca-e (21) 
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For the A(1232) case we replace the coupling constant NNtt (/) by the ANn 
coupling constant (/*) and the a and r 3 operators by the transition spin an isospin 
operators from 1/2 to 3/2 objects (with the normalization given by Eq. (§)). How- 
ever, the isoscalar part does not contribute in this case. Then, for the AN%° cou- 
pling we have: 

-tT =- e — CS ] -eTl (22) 
2 fi 

The vertices -fNNn , -fN*Nn° and 7 AiV7r of Eqs. flT|), (HJ) and @) are 
diagrammatically represented with the same Feynman graphs as the corresponding 
Kroll-Ruderman terms with charged pions in Fig. 3. 

The effective Lagrangians needed for the evaluation of the Feynman diagrams 
and the corresponding coupling constants are given in appendix A. The correspond- 
ing Feynman rules are given in Appendix B of ref. [pTj|. 

3 Results and discussions. 

In Figs. 5 to 10 we show the total cross sections for all the isospin channels, as well 
as the contribution to the total cross section of diagrams with nucleon, A(1232), 
AT* (1440), iV*(1520) and p(770) as intermediate states (in ref. [JTTJ] we also showed 
differential cross sections and invariant-mass distributions for the 7p — > n + n~p 
channel). We show results up to E n = 800 MeV, where the new experiments at 
Mainz concentrate, and we compare our results with the available experimental data 

g g, § §]. In the errorbars of the DAPHNE data [§ we have plotted only the 
statistical errors. The systematic errors are of the same order of magnitude. 

The isospin channels with one or two charged pions in the final state are shown 
in Figs. 5, 6, 8 and 9. We observe that the A(1232) terms (short-dashed lines) are 
dominant in these isospin channels. Essentially the A (1232) Kroll-Ruderman and 
A (1232) pion-pole terms (diagrams (i) and (k) of Fig. 3) are the more important 
terms. The non resonant terms (short-dash-dotted lines) are much smaller and they 
provide a small background which grows up moderately as a function of the energy. 
The iV*(1520) contribution (long-dashed lines) by itself is also small compared to 
the A(1232) one, but it is essential to reproduce the total cross section due to its 
interference with the A(1232) Kroll-Ruderman term as we already remarked in ref. 
[|TT|]. This interference occurs only between the A (1232) Kroll-Ruderman term and 
the s-wave part of the iV*(1520)A7r contribution. Since now part of the iV* — > An 
decay is due to <i-wave, the interference effects are smaller than in ref. [[n) where 
all the N* (1520) — > An decay was associated to s-wave. This has as a consequence 
a better agreement of theory with experiment than in ref. [ [TT| ] (see Fig. 5 here in 
comparison with Fig. 5 of ref. [[TT]]). The p-terms are negligible at these energies, 
but they show up clearly at higher energies (see ref. [|Tl|]). Contributions from other 
terms are still smaller. 

For the isospin channels with two neutral pions in the final state things are 
rather different (see Figs. 7 and 10). In these cases a lot of terms vanish due to 
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the fact that the photons cannot couple to neutral pions. In particular, there are no 
Kroll-Ruderman and pion-pole terms. Thus, the total cross section for these isospin 
channels is much smaller than the cross section for the other isospin channels. Fur- 
thermore, terms that in the other isospin channels are very small become important 
in these isospin channels. 

In Figs. 7 and 10 we can see that, except for the contribution from N-intermediate 
states (short-dash-dotted lines), which is very small, the other contributions are all 
of them relevant. Below 500 MeV we can see that the iV*(1440) (long-dash-dotted 
lines) dominates the reactions. At 500 MeV the A(1232) (short-dashed lines) and 
the iV*(1520) (long-dashed lines) start to grow up, and around 600 MeV all these 
diagrams have similar strength. At 600 MeV the iV*(1440) contribution starts to 
fall down, and the iV*(1520) and A(1232) dominate the reaction. The iV*(1520) 
contribution peaks at 720 MeV (it is responsible of the peak of the total cross 
sections) and from this energy on it falls down, while the A(1232) contribution 
continues growing up moderately. 

It is worth noting that, in these latter cases, the iV*(1520) contribution is im- 
portant by itself, and not by its interference with other terms as it happens in the 
isospin channels with one or two charged pions, and again it is essential to repro- 
duce the peak of the cross section around 700 MeV. To further stress the role of 
the iV*(1520) resonance, we plot in Fig. 11 the results for the 7p — ► ix G ix°p cross 
section using the two possible combinations of Eq. < ]T4| ) for the s- and d-wave 
couplings. We see that the curve (a) is clearly favoured by the data, precisely the 
same case which was favoured in the ^p — > n + n~p reaction. This provides further 
support for our choice of amplitudes in the A r *(1520) — ► Air decay. 

At this point, we should mention that in spite of the large iV*(1440) decay 
branching ratio into Nun (around 35%), and the fact that in the present reaction 
the iV*(1440) is on-shell at £ 7 ~ 600 MeV, the iV*(1440) contribution to the to- 
tal cross section is very small for some of the isospin channels. This is due to the 
relatively small iV*(1440) coupling to photons pD]], compared to the A(1232) or 
the iV*(1520) resonances (the branching ratio to N'j is around one order of mag- 
nitude smaller for the iV*(1440) than for the A(1232) or the TV* (1520)). However, 
it is important in the two-neutral pion channels in the low energy region (below 
E 1 ~ 650 MeV) (see long-dash-dotted line of Figs. 7 and 10). Note again, that for 
E 1 < 600 MeV the iV*(1440) is the dominant contribution, and although there are 
some discrepancies with de DAPHNE data in this region, our model agrees well 
with the preliminary data of the TAPS collaboration in the same region. 

For the total cross sections, in Figs. 5 to 10 (solid lines) we can see a remarkable 
agreement with the experiment in some isospin channels, and some discrepancies 
in others. 

For the ^p — > n + ix~p channel (Fig. 5) we reproduce quite well the total cross 
section up to = 800MeV^ (we also reproduce invariant mass distributions even 
at higher energies [[TT]]). 

For the ^p — > ir + 7r°n channel (Fig. 6) we have found an important discrepancy 
between our calculations and the experimental data. It is very easy to understand 
qualitatively our results for the 7p — > n + n°n cross section from isospin coeffi- 
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cients. As we have already said, the photoproduction of one or two charged pions 
is dominated by the A (1232) Kroll-Ruderman and A(1232) pion-pole terms. In the 
7P — ► n + n~p channel, these terms with a A ++ in the intermediate state have a fac- 
tor 1 in the amplitude from isospin (T A ), while in the 7p — > n + n°p case with a A 
in the intermediate state, this factor is y/2/3. Then we have that the cross section 
for the 7p — > n + n°n channel is around 4.5 times smaller than the 7p — > n + n~p one 
(this relation is not exact because there are other terms with different isospin coef- 
ficients that also contribute). As we can see comparing Figs. 5 and 6 this relation is 
approximately satisfied. 

In order to understand these discrepancies, we have tried to look for additional 
contributions for this channel. For instance, we have calculated the contribution to 
the total cross section of new diagrams with p-intermediate state decaying into WW 
where the photon is coupled to the NNp vertex in one case and to the intermediate p 
meson line in another case. Note that these additional diagrams only contribute for 
the 7£> — > n + 7r°n and yn — > 7r~7r°p channels, because for the others channels, the 
intermediate p-meson is neutral, and does not couple to photons. However, we have 
found that these contributions are very small, and they do not modify our results. 

On the other hand, we have also investigated the contribution to the total cross 
section of the off-shell part of the A interactions with photons and pions (see ref. 
[11H]). For this purpose, we have evaluated diagram {%) of Fig. 3, including the 
off-shell part of the A -interaction. Then, from the off-shell parameters of the A in- 
teraction, we get an additional term which is very similar to the one obtained from 
diagram {%) of Fig. 3 without the off-shell parameters, but it has a multiplicative 
factor (y/sA — Wa)/2wia (y/sK is the A invariant mass, and m& is the A mass). 
This factor is zero for a A on-shell, and in our range of energies, it changes between 
—0.05 and +0.07, and it strongly peaks around zero (see Fig. 6 of ref. [TTTp). This 
is so, because the pion of the 'jANit vertex in Fig. 3 (i) picks up the necessary 
momentum to make the A closest to on-shell. When we perform the actual calcula- 
tions the modifications from the consideration of these off- shell corrections are of 
the order of 1% for a broad range of the off- shell Z parameter (we have checked 
it for several values of Z between —1/2 and 1/2). Then, we conclude that the off- 
shell contribution of the A(1232) resonance to the total cross section is very small, 
and we can neglect it. 

For the 7p — > 7r°7r°p channel (Fig. 7) we can see a good agreement with the 
experimental data above 650 MeV, but our results are below the experimental ones 
for E 7 < 650 MeV. We should also mention that there are already preliminary 
results for this channel from the TAPS collaboration [|9|], which for E 1 < 600 MeV 
are below the DAPHNE data, and our model agrees very well with these preliminary 
results in this energy region. Nevertheless, we should wait for their final results. 

For the •yn — > n + n~n we can see in Fig. 8 that we approximately reproduce the 
experimental data below 650 MeV, but above this energy our model overestimates 
the experimental data. 

For the yn — > 7r~7r°p isospin channel (Fig. 9), the theory is a bit below the ex- 
perimental results [[3J, 0fl, but the trend of the energy dependence is well reproduced. 

However, in connection with the data of refs. [0, |]] for the yn — > 7r + 7r~n 
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and jn — ► 7r~7r°]9 reaction, we should point out that there is some disagreement 
between the results of the same experiments for the 7p — > n + n~p channel and the 
experimental results of ref. [[[], [8|] (see Fig. 5). One should then be cautious and wait 
for new measurements in these channels before extracting any conclusion from the 
present discrepancies of our model with these data. 

In Fig. 10 we show the results of our model for the yn — > n n°n channel for 
which there are not yet any experimental data. 

Due to the interest raised on the two pion photoproduction at threshold [T3], 



14, 15, 16], we have also calculated the total cross section for the yp — ► n + n p 



and 7p — > ir n°p channels at energies near threshold. In Figs. 12 and 13 we show 
our results for the total cross section for these channels, compared to the results of 



ref. Q13J, |14J (in ref. D13Q only the amplitude at threshold is given, not the total cross 
section. Then, we have taken that amplitude as a constant, and we have integrated 
it over the three body phase space). 

In Fig. 12 we see the results for the yp — > n + n~p channel. There appear to be 
some discrepancies between the results of refs. JT3| ] and [ |14| ] which are mostly due 
to the fact that we use the constant amplitude at threshold of ref. [O. Indeed, if the 
same is done for the work of ref. JT4"| ] the results are very similar (see Fig. 9 of ref. 
[|14|]). This also means that energy dependent terms are important even very close 
to threshold. 

Our results are above the results of refs. [ p"3"| , [14| ] . We should note that in addition 
to the Chiral terms which contribute at threshold, and which we have in terms of 
our effective Lagrangians, there is one important contribution at threshold which 
is missed in the two other approaches. This corresponds to the diagram of Fig 3. 
(q) and the crossed one which is not drown there but is actually calculated in this 
channel. They correspond to iV* (1440) excitation with decay into nn in s-wave and 
N. These terms were also shown to be important at threshold in the nN — > nnN 
reaction [[18]. Its relevance in that reaction has been stressed in a recent paper [ j32[ ] 
where the existence of an additional energy dependent term of unknown strength is 
pointed out. 

The different threshold behaviour of the approaches of refs. [ |T3| , [T4] ] and the 
present work, apart from some differences in the treatment of the A terms, is mostly 
due to the fact that we include the iV*(1440) resonance. In addition there is a small 
contribution from the iV*(1520) a few MeV above threshold, and there are also 
other terms which vanish at threshold but contribute close to it. 

We should also mention that in the region of interest to us, 450 MeV < E 1 < 
800 MeV, the cross section is dominated by resonance terms which are not pro- 
vided by Chiral Perturbation Theory. 

In Fig. 13 we show the threshold results of the •yp — > n°n°p channel. We 
observe large discrepancies between the results of refs. [ P"3"| ] and [ |T4"| ] which are 
commented in refs. [ [T3i Among other differences, the results of ref. [[14]] 
include loop corrections which are relatively important in this channel close to 
threshold [jnfl. Once again our results, which also miss the loop corrections, are 
bigger than in the two approaches of refs. [ [TS] , [Bj] for the same reasons discussed 
above. Our results at threshold also include the contribution of the crossed term of 
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the iV*(1440) diagram of Fig. 3(q), which is negligible at higher energies. 

Furthermore, we should also stress that, like in the 7p — > n + n~p channel, the 
terms which vanish at threshold dominate the reaction at large energies. 

One limitation of our model is that we neglect unitarity. Unitarity for three par- 
ticles in the final state is a nontrivial problem [J33|], and is not the purpose of this 
paper to deal with it. However, for the dominant channel (7 A — > 7rA[7rA]), due 
to the fact that the A-width is implemented in the A-propagators, our model ap- 
proximately satisfies unitarity. In order to estimate the errors that unitarity could 
introduce, we follow the procedure of Ols son [ [34"| ] to unitarize this channel by mul- 
tiplying the A-terms by a phase and requesting that the resulting amplitude, after 
adding the background, has the phase of the tcN amplitude. The angle <\> needed 
for this phase, e 1 ^, is of the order of 10°. We have checked, for the 7p — > n + n~p 
channel, that implementing this phase in the amplitude changes the results at the 
level of 3% at energies below E 1 = 800 MeV. Even increasing the angle <p to 
20° the changes are of the same order of magnitude. This result is due to the fact 
that, except for the already commented interference between the iV*(1520) and the 
A Kroll-Ruderman term, the interference between the different terms of the ampli- 
tude is small. 

This is only a rough approach to such a difficult problem, however it gives hints 
that the unitary corrections might be small in that energy range. This conclusion 
is also supported by another approximate scheme to unitarize the Born amplitudes 
used in [jllj. An absorptive correction factor is used there which does not modify 
the cross section below = 800 MeV, although it produces a substantial reduc- 
tion at much higher energies. Nevertheless a rigourous treatment of this problem 
would be welcome. 



4 Conclusions. 

We have constructed a model for the 7 A — > nnN reaction extending the model 



of ref. JTip for the jp — > ix+ix~p reaction to the rest of the isospin channels, 
including some improvements in the A(1232), iV*(1440) and iV*(1520) couplings. 
Our model reproduces quite well the experimental results of refs. below 
E 1 = 800 MeV for the jp — > n + ix~p and jp — > 7r°7c°p isospin channels, but we 
have found some important discrepancies for the ^p — > n + n°n channel. For the 
isospin channels on the neutron we have found also some discrepancies with the 
old data of refs. [0, ^p, but as we have already mentioned, these experiments should 
be improved. 

As we already observed in ref. [ [TT| ] the reaction is dominated by the A(1232) 
Kroll-Ruderman and A(1232) pion-pole terms, but we get an appreciable contribu- 
tion from other terms. In particular we have shown the crucial importance of the 
A'"* (1520) terms. In the isospin channels with one or two charged pions in the final 
state, the A* (1520) contribution is small compared with the A(1232) terms, but it 
shows up clearly in the total cross section due to its interference with the A(1232) 
Kroll-Ruderman term. For the cases with two neutral pions in the final state the 
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dominant terms in the other isospin channels vanish, and the contribution of the 
iV*(1520) terms is bigger than the contribution of the A(1232). At energies be- 
low E 1 = 600 MeV the iV*(1440) plays an important role in the two neutral pion 
channels. 

One of the interesting findings was the possibility to extract information on the 
coupling constants for iV*(1520) decay into Air. The interference of the iV*(1520) 
excitation term, followed by A7r decay, with the dominant A(1232) Kroll-Ruderman 
term in the yp — > n + n~p channel allowed us to determine the sign of the couplings 
for the s-wave and d-wave in the iV*(1520) — > Air decay, when the absolute val- 
ues given by the branching ratios, and the constraint A s /Ad > obtained form the 
ttN — > titxN reaction were used. We also observed that the same choice of ampli- 
tudes gave rise to a fair agreement with the yp — > 7c°ir°p channel while the other 
possible choice gave rise to a cross section in sheer disagreement with the data. 

While in the •yn channels there are reasons to expect that the data will change 
when new experiments are done, the discrepancy found with the recent data pSp 
of the 7p — > Ti + 7i n channels is puzzling. For the moment we cannot envisage 
a solution to this problem. Alternative experiments would be in any case most 
welcome. 



Acknowledgements: Discussions with M.J. Vicente- Vacas are warmly acknowl- 
edged. We would also like to thank N. d'Hose, G. Audit, G. Tamas and H. Stroher 
for multiple discussion around their preliminary data. This work has been par- 
tially supported by CICYT contract number AEN 93-1205. One of us J. A. Gomez 
Tejedor wishes to acknowledge financial support from the Institucio Valenciana 
d'Estudis i Investigacio. 



14 



LAGRANGIANS . 



Appendix A 



Cnn. = — V-y^dJ ■ r V (A.l) 
fx 



Ai— - - 

Cnn™ = -4vr— ^0 • 0^ (A.2) 
H 



f* 



Can, = —VlSiidi^T^N + h.c. (A3) 



Ua, = — W A )^tt A + h.c. (A.4) 



— f -> 

jCn*Nk = —V N *°i (dif) ■ t ^ N + h.c. (A.5) 



Cn*n™ = -CV N *<f>-<f>V N + h.c. (A.6) 



C N *An = -9^-tf A Si(di(i> x )T x * N . + h.c. (A.l) 



Cn>*a. = Mn» f - 9j ^S\d l SjdA A T A ^ A + h.c. (A.8) 
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-eVNi-fAp - ^a^d u A^ N (A.9) 



s(<M"0- - <M M 0+K (A.10) 



-iq- ^fflsAj -ftf (A. 11) 



-iq— *^^A0 A ^ A *jv + /i.e. (A.12) 
/i 



f + 

-iq-q#V*M*0 • "r*N + /i.e. (A.13) 
/i 



^2*^ e . jfc5 -t (9 . Afe) Wjv + (A14) 



-^-IfNa^duA^N* + h.c. (A. 15) 



[g-,S • A - ig a (a x • + /i.e. (A.16) 



{^ J v,7 M 0!f ) - %^^0(f)) • (A.17) 



2m 
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C p ^ = -fjf -(4>xdy) (A.18) 



C ANp = -^ P —^ikSi (drftf^ T x m + h.c. (A.19) 



£>n'*N(P = -g P ^ N Si4 P) ^N'* + h.c. (A.20) 



£ p7TJ = imfthi daApfdLtfP (A.21) 
A 4 



0^- 7 = e/ p <^ tt+A^- + 0-^0+) (A.22) 



£™r 7 = ^-e^e abc A, d v (f> a d a( j> b d fi </> e (A.23) 



L V11 = El e ^<f F ^ Fap (A.24) 

Ffiu = dfj.Au — d v A^ 



In these expressions, 0, \I>a> ^tv*> ^jv*. 0|f an d A M stand for the nucleon, 
pion, A(1232), iV*(1440), JV*(1520), p(770) and photon fields, respectively; m 
and /U are the nucleon and pion masses. 

For the A A7 coupling we write directly the vertex contribution to the Feynman 
rules [ f35[ ] in analogy to the NN^y case: 

- i5H^ = »( fr + p/) e A + « x fc)l • e (A.25) 

I 2m^ 0111 I 



17 



where Sa is the ordinary spin matrices for a spin 3/2 object, p and pf stand for the 
initial and final delta momentum respectively, and k for the photon momentum; 
and e& are the magnetic moment and charge of the delta; is the delta mass. 
The coupling constants are listed below. 



Coupling constants: 

f = l /* = 2.13 
A x = 0.0075 

f A = 0.802 /ajv 7 = 0.116 

/ = 0.477 C = -2.29/i- 1 

9an*k = 2.07 e = 0.3027 

/at'-Att = 0.911 ^Att = -0.552 



1.79 for proton ~ N _ j 0.0173 for proton 

—1.91 for neutron 7 { —0.0112 for neutron 



„ N J 0.108 for proton „ N j —0.049 for proton 

^ 7 1 —0.129 for neutron ® a 1 0.0073 for neutron 

G^NNp = 2-9 G'nnp = 18.15 

/ P = 6.14 C p = 2 

# p = 0.591 # p7r7 = 0.0378 

F 37r = 0.0259/i- 3 P = 0.0035/x- 1 

t? = V a*a++ = 1-62 Aip 

The magnetic moment of the A(1232), a^a* can be calculated in the quark model 
[[361], w i m m e result //a//-^ = e A/e. We shall use this result, except for the A ++ 



where we shall use the experimental value /xa++ = (1-62 ± 0.18) pL p based on the 
rip bremsstrahlung (ir + p — > n + p'j) [|37|]. 
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Appendix B 



TOTAL CROSS SECTION. 

The cross section for the 7^ — > irirN reaction is given by 



m 

a 



1 r d 3 p4 r d 3 p 5 r 3 m 
VkY I ~2uu J ~2um J P2 E~ 9 



A 1 /2( s? o,m 2 ) (2tt) 

5 4 (A , + pi _ p2 _ pi _ p^Y, \T\ 2 (B.l) 



Si Sf 

m 2 1 



doj^duiid cos 6*5(i045 



~~ A 1 /2( S; o,m 2 ) 4(2tt) 4 



Where k = (u,k),p 1 = (E 1 ,p 1 ),p 2 = (E 2 ,p 2 ),Pi = (uj 4 ,p 4 ),p 5 = (u 5 ,p 5 ) 
are the momenta of the photon, incident proton, outgoing proton and the outgoing 
pions respectively. In (B.2) 45 , #45 are the azimuthal and polar angles of p^ with 
respect to p$ and # 5 is the angle of p§ with the z direction defined by the incident 
photon momentum k. T is the invariant matrix element for the reaction. We repro- 
duce this formula here since in Eq. (2) of ref. [ |TT| ] there was a missprint and d6$ 
should have been dcos9 5 , as in Eq. (B.2) above. 
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Figure captions: 

Fig. 1: Classification of the Feynman diagrams into one point, two point and 
three point diagrams. Continuous straight lines: baryons. Dashed lines: pions. 
Wavy lines: photons and p-mesons (marked explicitly). 

Fig. 2: Total cross section for the 7p — > it + it~p reaction, (a) fN'*Aw = 0.911, 
gN>*Air = -0.552; (b) /a^att = -0.911, g N >* A K = 0.552. 

Fig. 3: Feynman diagrams for the 7p — > n + n°n, 7p — > n°n°p, jn — > n + n~n, 
7n — > TT~n°p and — > 7T 7T°n reactions. 

Fig. 4: Feynman diagrams for the *yN — > n°N amplitude. 

Fig. 5: Total cross section for the 7p — > n + n~p reaction. Continuous line: 
total cross section. Short-dashed line: contribution of A(1232)-intermediate states. 
Long-dashed line: contribution of iV*(1520)-intermediate state. Short-dash-dotted 
line: contribution of iV-intermediate states. Long-dash-dotted line: contribution of 
p-intermediate states. Short-dash-long-dashed line: rest of the diagrams. Experi- 
mental data from refs. [jl], 0, |, |, §. 

Fig. 6: Total cross section for the 7p — > n + n°n reaction. Continuous line: total 
cross section. Short-dashed line: contribution of A(1232)-intermediate states (dia- 
grams (i)-(o) of Fig. 3). Long-dashed line: contribution of iV*(1520)-intermediate 
state (diagram (p) of Fig. 3). Short-dash-dotted line: contribution of N- and p- 
intermediate states (diagrams (a)-(h) of Fig. 3). Long-dash-dotted line: contribu- 
tion of N*( 1440) -intermediate states (diagrams (q)-(t) of Fig. 3). Experimental 
data from ref . pSp . 

Fig. 7: Total cross section for the 7p — > ir n°p reaction. Continuous line: total 
cross section. Short-dashed line: contribution of A(1232)-intermediate states (dia- 
grams (i)-(o) of Fig. 3). Long-dashed line: contribution of iV*(1520)-intermediate 
state (diagram (p) of Fig. 3). Short-dash-dotted line: contribution of iV-intermediate 
states (diagrams (a)-(g) of Fig. 3). Long-dash-dotted line: contribution of N*( 1440)- 
intermediate states (diagrams (q)-(t) of Fig. 3). Experimental data from ref. [|8p. 

Fig. 8: Total cross section for the — > n + n~n reaction. Continuous line: total 
cross section. Short-dashed line: contribution of A(1232)-intermediate states (dia- 
grams (i)-(o) of Fig. 3). Long-dashed line: contribution of iV*(1520)-intermediate 
state (diagram (p) of Fig. 3). Short-dash-dotted line: contribution of N- and p- 
intermediate states (diagrams (a)-(h) of Fig. 3). Long-dash-dotted line: contribu- 
tion of iV*(1440)-intermediate states (diagrams (q)-(t) of Fig. 3). Experimental 
data from refs. [0,0. 
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Fig. 9: Total cross section for the — ► n~TT°p reaction. Continuous line: total 
cross section. Short-dashed line: contribution of A(1232)-intermediate states (dia- 
grams (i)-(o) of Fig. 3). Long-dashed line: contribution of N*( 1520) -intermediate 
state (diagram (p) of Fig. 3). Short-dash-dotted line: contribution of N- and p- 
intermediate states (diagrams (a)-(h) of Fig. 3). Long-dash-dotted line: contribu- 
tion of N*( 1440) -intermediate states (diagrams (q)-(t) of Fig. 3). Experimental 
data from refs. [|,|p. 

Fig. 10: Total cross section for the — > 7r°7r°n reaction. Continuous line: total 
cross section. Short-dashed line: contribution of A(1232)-intermediate states (dia- 
grams (i)-(o) of Fig. 3). Long-dashed line: contribution of A^*(1520)-intermediate 
state (diagram (p) of Fig. 3). Short-dash-dotted line: contribution of iV-intermediate 
states (diagrams (a)-(g) of Fig. 3). Long-dash-dotted line: contribution ofiV* (1440)- 
intermediate states (diagrams (q) and (t) of Fig. 3). 

Fig. 11: Total cross section for the 7p — ► ir n°p reaction, (a) fa'* An = 0.911, 
9n'*Att = -0.552; (b) f N »± K = -0.911, g N >*An = 0.552. 

Fig. 12: Total cross section for the ^p — > n + n~p reaction at energies close to 
threshold. Continuous line: results of our model. Long-dashed line: results of M. 
Benmerrouche and E. Tomusiak JT3|]. Short-dashed line: results of V. Bernard et al. 

[0- 

Fig. 13: Total cross section for the ^p — > 7r°7r°p reaction at energies close to 
threshold. Continuous line: results of our model. Long-dashed line: results of M. 
Benmerrouche and E. Tomusiak [jT3|]. Short-dashed line: results of V. Bernard et al. 
& 
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